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ABSTRACT: An electron density functional approach for the calculation of the 
nuclear multipole moments is presented. The electronic matrix elements entering 
the experimentally observed hyperfine electron-nucleus interaction constants in 
atoms are expressed in terms of the electron density functionals of the charge 
or spin distributions. In principle, the construction of the charge or spin dis- 
tribution density functionals can be obtained by means of every relativistic or 
non-relativistic quantum-mechanical or DFT method. The electronic matrix el- 
ements for all the electronic operators of the existing Hamiltonians of hyperfine 
electron-nucleus interactions are considered. 
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1 INTRODUCTION 

The hyperfine structure in the spectra of the many-electron systems is a result 
of the interaction of the electronic and nuclear multipole (electric and magnetic) 
moments. This interaction leads to the superfine splitting of the electronic spec- 
trum, the rotational spectrum of molecules and the electron paramagnetic reso- 
nance lines as well [1, 2]. The energy of the hyperfine electron-nucleus interaction 
is determined by the finite nuclear size, the nuclear spin orientation and the nu- 
clear magnetic moments space distribution. This energy is much smaller than 
the energy of the interactions in an electronic system and must be considered 
separately and independently for each fine structure component. 

The transitions between neighbouring levels of the fine and hyperfine struc- 
tures can be measured by means of highly sensitive radioscopic methods, for ex- 
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ample, microwave absorption, electron-spin resonance (ESR), nuclear magnetic 
resonance (NMR), atomic and molecular beams. 

These experimental methods [1, 2] can give us information about many prop- 
erties and effects of nuclei connected with nuclear magnetism, for example: the 
mechanism and time of spin relaxation, the width of the nuclear spin resonance 
lines and the indirect interaction between nuclear spins. Using the methods of 
NMR we can experimentally determine the magnetic nuclear moments. There is 
no any experimental method for direct measurement of the electric nuclear mo- 
ments and especially the quadrupole momentum. They can be determined with 
great accuracy from the hyperfine splitting in electronic and molecular spec- 
tra [3] -[5]. For some heavy elements there are some good muonic values of the 
quadrupole momentum. 

In the analysis of the multipole electron-nuclear interactions the Hamiltonian 
is represented as a sum of products of the electronic and nuclear operators. These 
operators can be written down explicitly and can be obtained by using the per- 
turbation approach by expressing the interaction Hamiltonian in series in powers 
of the ratio of nuclear and electronic coordinates [3, 6, 7], as well as by a direct 
introduction of all the multipole interactions [8, 9]. 

The operators of the magnetic electron-nuclear multipole interactions can be 
determined also by expanding the Hamiltonian in series in powers of the constant 
of the fine structure which leads to some relativistic corrections in the main 
Hamiltonian [10]- [12]. 

Usually the experimental hyperfine multipole interaction constants are repre- 
sented as a product of the matrix elements of the corresponding electronic and 
nuclear operators [1, 3, 6, 7]. The diagonal matrix elements of the nuclear opera- 
tors correspond to the multipole moments (magnetic dipole, electric quadrupole, 
magnetic octopole nuclear moments) [1, 6, 13]. In addition to the diagonal elec- 
tronic matrix elements, the non-diagonal matrix elements must be considered as 
well, because of some important second order effects. The non-diagonal matrix 
elements are also important in the analysis of the magnetic transitions induced 
by the external magnetic fields [3]- [5]. Using the theoretical interpretation of the 
multipole electric and magnetic moments through the experimental values of the 
hyperfine interaction constants and the calculated multipole electronic operators 
matrix elements, one can obtain the values of the electric and magnetic nuclear 
moments. The calculation of the matrix elements of the multipole magnetic 
operators is essential in determining the relative contribution of the magnetic 
interactions and gives us the possibility to estimate the role of the quadrupole 
interaction in a hyperfine structure. Different transitions in the resonance spectra 
make it possible to extract the quadrupole interaction constant by means of the 
analysis of the magnetic hyperfine sub-term or by means of the analysis of the 
deviation from the interval rules produced by a magnetic dipole electron-nuclear 
interaction. The correct description of the hyperfine constants needs the inclusion 
of the relativistic corrections [14]. 
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In this paper we present the electronic matrix elements for all electronic op- 
erators and express them in terms of the density functional of charge or spin 
distributions. Thus for including the relativistic effects one must take the corre- 
sponding density functional of relativistic corrections, for instance, the functional 
given in paper [15, 16]. 



2 ELECTRON-NUCLEAR OPERATORS OF THE 
HYPERFINE INTERACTIONS 

2.1 Multipole Interaction Operators 

The operators of electron-nuclear multipole interactions can be written down 
[6, 7] 

h = EH fe = EQ" • = EE(- 1 W- C 1 ) 

k k k q 

where H fc is the Hamiltonians of the electric (fc-even) and of the magnetic (/c-odd) 
interactions of the nuclei with the electronic shell. Here Q fc and F k are irreducible 
tensor operators of k rank with components Q k and F k , q = (—k, —k + 1, k — 
l,k). These operators act in the spaces of the total nuclear moment I and of the 
total electronic moment J correspondingly. 

The electronic operators F k in the interaction Hamiltonian expressed in series 
in powers of space coordinates are represented in the form 

, u [ -E,rr (fc+1) C fc (^), {k - even) 

F k = I , (2) 

, Ei^eV i rr (fc+1) C fc (0 i ,^)(|L(i)-2S(i)), (A; -odd) 

where \i e is the Bohr's magneton, L(i) and S(i) are the orbital and spin angular 
momenta of the electron with number i, = (r i5 ifii) is its radius-vector and 
C k (6i,(fi) are the normalized spherical functions. 

The diagonal matrix elements of the Qq component of the nuclear operator 
Q fc in the state /, Mi = I, is the 2 h total nuclear electric (even k) or the nuclear 
magnetic (odd k) moments 

(II\Q k o\U) = ^=(7/*^ </||Q fc ||/>= 

M k (k - odd) 
Q k (k - even) 

where C\ \ I I are the Clebsch-Gordon coefficients and < J||Q fc ||J > is the re- 
duced matrix element of the operator Q fe . The matrix elements with k — 1, 2, 3 



3 



correspond to the nuclear moments [6, 13] 



m 1 //. a. |, m 3 = -ft, 

where \i is the magnetic dipole, Q is the electric quadrupole, and ft is the magnetic 
octopole nuclear moments. 

Let us consider the case when the external field can be neglected. Then in 
the first approximation of the perturbation theory the energy of the hyperfine 
interactions is determined by the matrix elements of the Hamiltonian H Eq.(l) 
which are diagonal with respect to J, I and total moment F = I + J 

E ( j ] (F) = < (3IaJFM F \U\(3IaJFM F > . (4) 

Here M F is the value of the z- projection of the total atomic moment, a and 
(5 are the additional quantum numbers of the electronic and nuclear systems 
correspondingly. 

Using Racah formalism [17] we can rewrite the expression in the form 

E { j\f) = ^(-l) /+J - F H/(//JJ; A;F) < aJ\\F k \\aJ >< f3I\\Q k \\f3I > 

k 

= *£A k M(IJ;F;h), (5) 

k 

where 

A k = < aJ J\F k \aJ J >< (3II\Q k \(3II >; 

M(/J; F; k) = (-1) I+J ~ F — W(IIJJ; kF), 

dk 

_ \ r<j k j r<i k i 

a k — r — J °/ Ii 

V /(2J+1)(2/+1) 

and W(IIJJ] kF) are the Racah coefficients. 

The first three hyperfine interaction constants are connected with the 

experimental values of the dipole magnetic (A), the quadrupole electric (B) and 
the octopole magnetic (C) interaction constants as follow 

Ax = IJA, A 2 = 1/45, A 3 = C 

Calculating the electronic matrix elements in Eq. (5) and using the experimental 
values for Ej~\f) and Ak we may obtain the nuclear moments /i, Q and ft. 
The second order perturbation term 

F (2W ^ | < (3 la J FMp \~H\f3Ia J'FMp > | 2 

fi j (^) = L p — p > w 
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where Ej is the energy of the fine structure level J, becomes important in the 
case when the distance between hyperfine structure levels is of the same order as 
the energy intervals of the fine structure. 

In this case the constant corresponding to Ej '(F) in Eq.(5) must be added, 
as a correction to the hyperfine interaction constant Ak . 

2.2 Breit-Pauli Electron-Nucleus Operators 

The Hamiltonian H en of the electron-nuclear interactions obtained by means 
of the expansion of the electronic system Hamiltonian in series in powers of 
the fine structure constant includes three operators of the effective Breit-Pauli 
Hamiltonian [10]- [12] 

H en = + Hf^ + Hf/ nt . (7) 
The electron-nuclear dipole-dipole interaction has the following form 

Ug = g aga P Y: {3[S(<) • r,][I • r 4 ] - r?S(i) • i} rf. (8) 

i 

Here g$ is the ^-factor of the free electron, a = l/2c is a half of the fine structure 
constant (in atomic units), g is the ^-factor of the nucleus and a p = l/2m p c. The 
magnetic interaction between the electron angular momentum and the nuclear 
moment has the form 

H£ = 2a<7ap5>f S L(0-I. (9) 

i 

The operator of the electron-nuclear contact interaction is 

H£* = ^goaga p J2S(r t )S(t) • I. (10) 

6 i 

It is straitforward to point out that the direct inclusion of the hyperfine interac- 
tions of the nuclear momenta and the electronic system leads to the same results 
as the Hamiltonian H en in Eq.(7) [8, 9]. 

3 MATRIX ELEMENTS OF THE ELECTRON 
OPERATORS 

In order to make the matrix elements of the electronic operators suitable for cal- 
culations we have to recouple them from the LMlSMsJMj coupling scheme to 
the LM L SM S coupling scheme. Since we consider the hyperfine splitting sep- 
arately for each level of the fine structure, applying Racah techniques [17] we 
express the matrix elements < aLSJ\F k \aLSJ' > and < aLSJ\tl en \aLSJ' > 
directly through the diagonal matrix elements of F fe and H en in the LMlSMs 
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coupling scheme. Thus, along with the diagonal in J matrix elements determin- 
ing the A)~ constants,we have to calculate the non-diagonal in J matrix elements 

f<y\ 

determining A k values. It must be pointed out that the non-diagonal in J ma- 
trix elements are very important in the estimation of the transition probabilities 
between the hyperfine levels [2, 4]. 

3.1 Electric Multipole Operators 

The symmetrical tensor operators (F el ) k Eq.(2) with the components 
(F e ')/c _ ^ i — r^ k+1 ^Cq (6i, ifi) act only on the space variables, assuming the 
coupling scheme corresponding to J = L + S. Then following the Racah technique 
[17], and using the diagonal matrix elements (F el ) k in the LM L SM S coupling 
scheme we may express the electronic matrix elements of Eq.(6)as 

< aLSJ\F k \aLSJ' > = (-l) s ' L - J+k+1 ^ [2J + 1)(2 J' + 1) 

x W(LLJJ';kS) <aLS\Y l -r7 {k+1) C%(6 i ,<p i )\aLS> . (11) 

i 

Taking into account that the operators r~^ k+1 ^C k (6i,ipi) are operators of mul- 
tiplication the diagonal matrix elements can be expressed by means of the lo- 
cal density functional of charge distribution p{KK\v) = p{K\v) in state K = 
\aLM L SM s > in the form 

< aLS\ J2r- {k+1) C k (6 t ,^)\aLS > = j r-^C k {9^)p{K\v)dv. (12) 

i 

Up to now there is no evidence of the existence of the electric electron-nucleus 
interactions with k > 2. The quadrupole interactions (k = 2) can be observed 
only in the case of degenerate nuclear states and nuclear spin I > 1 [1]. In our 
approach in the case of the quadrupole interaction the electronic system may 
posses an arbitrary total moment. As a matter of fact in the case under consid- 
eration the density charge distribution does not depend on the electronic system 
degeneration. Calculation of the matrix elements of the electronic operator of 
the quadrupole interaction is equivalent [??] to the calculation the matrix ele- 
ments of the gradient of the electrostatic field created from the electrons placed 
in the centre of the nucleus [l],[3]-[5], [7]. 

3.2 Magnetic Multipole Operators 

For each value of k the operators Eq.(2) of the magnetic multipole interactions 
are the sum of two operators. The first term acts only on the space variables of 
the electronic system 

(K 9 ) k = 2 ~y E[v,(^ (fc+1) c^, n))] ■ L(i). (is) 
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In terms of contravariant and covariant cyclic components this operator may be 
written as 

(K 9 ) k [v^ k+1) c\e^ riTLlii). (14) 

K i m=-l 

The contravariant component of the vector in the middle brackets, corresponding 
to the component C k (9i,<fi) of the normalized spherical function C k {6 i) ip i ), can 
be written as 



[V,(r- (fe+1) Cj(^,^))] m = ^y^[^(r- {k+1) Y kq (9 h ^)T (15) 

= ^(k + l)(2k + 3)r^ k+2) C k ^ I JCj^, <pi), 

where Y^ g hl (6 ) , ip) are spherical vectors [18], and Cj^ ^ J are Clebsch-Gordon 
coefficients. The scalar product in the right side of the expression Eq.(14), cor- 
responding to the component C k (9i,(fi) takes the form 

T,[^(n (k+1) c k q (e t , Vi ))] m Ll n ( i )= (16) 

m 

J(k + l)(2k + 3)r-( k+2) [C k+1 (i) x L(0]}, 

where 

[C fe+1 «XL«]J = ECj^mSC^i^i^mW- 

m 

Hence, the operator (F™ 9 ) fc is an irreducible tensor operator olk — th rank. Since 
the operators (F™ g ) k act only on the space part of the electronic system, for the 
matrix elements of the q component (F™ 3 )^, we obtain 

< aLSJ\(F™ q ) k q \aLSJ' > = (_i)S-w+k+i 

x ^(2 J + 1)(2 J' + l)(k + l)(2k + 3)W(LLJJ'; kS) £ C£» I \ 
x < aLS\ J2r7 (k+2) [C k+1 (i) x L(i)]J|a;LS > . (17) 

i 

We can represent the sum in Eq.(17) as 

= E^K.m! / rr (fc+2) c;^(^, ^1)^(1)^1^;^)^, (is) 
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where p(K\r;r') is the non-local charge distribution density functional in 

K = \aLM L SM s > state. In Eq.(18) the density functional does not depend on 

the summing indices m and q. 

The second operator of the magnetic multipole interactions acts on the space 
coordinates, as well as on the spin components of the electronic system 

(FD fc = 2/i e E Vitrr^.w)] • S(<). (19) 

i 

Now we obtain 

(FD fc = 2^(k + l)(2k + 3)J2ri {k+2) [C k+ \ l ) x S 1 (^] fc , (20) 

i 

The q component of the tensor product in Eq.(20) has the form 

\c k+i (i) x s(oi; = m ; c^e^sM), (2i) 

m 

where the spin moment acts on the spin subsystem. Since both operators C k+1 (i) 
and S(i) act on two different subsystems they commute , and the matrix elements 
of the components (F™ 9 ) k can be expressed as 



x 




< aLS J\(F™ g ) k \aLS J' > = 2p e y(k + l)(2k + l)(2k + 3)(2J+ 1)(2J' + 1) 

> ^yfc+l 1 fe 
/ j q—m m q 
m 

x < aLS\ tt^CWfariSMlaLS > . (22) 

i 

Applying the equations Eqs.(61)-(64) of the Appendix to the right hand matrix 
elements in Eq.(22) we have 

< aLM L SMs\Y^ r \ ik+2)c q-nMi^ ^i)Sl n (i)\aLM L S M' s >= 

i 

J r-^^^^Og^lr^ri)^ = 

r'j=ri 

< SMslS'jSM's > J r^ k+2) C k +U$i, ¥>i) (23) 

where D s (ri) = q(K K\ri, ri)l/S is the normalized spin distribution density 
functional, K' =| aLM L SM' s >, and K =| aLM L SS >. 

The application of the Wigner-Eckart theorem for the matrix elements < 
SM s \Si\SM' s > yields 

< SMslS^SM'g >= ^^^fMs^h- ( 24 ) 
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Finally the contribution of the component (F™ 9 )q of the magnetic operator F™ 9 
for the hyperfine splitting of any spin multiplet term under consideration (SMs) 
has the form 



< aLSJ\(F™ 9 ) k n \aLSJ' > = V6fi e - 1 = 
x yj(k + l)(2fc + 1)(2A; + 3)(2J+ 1)(2J' + 1) 

r l l fc+i I 

v J c q i (vrw i fc^ 51 1 5 

x S D D 1 ( Z^ U q-m m q U M s mM' 



x J r-^C^^Ds^dv,. (25) 



3.3 Electron-Nucleus Dipole Operator 

The magnetic dipole coupling between electron and nuclear spins Eq.(8) can be 
written down in a form in which the electronic and nuclear variables are separated 

H£ = 9oaga p £ S(i) • [3^ - ^] • I. (26) 

Here n 2 (i) is the unit tensor formed from the unit vector n(i) = r^/V;. 
The sum of the one-electronic operators in Eq.(26) can be written as 

HfJ = 0oa0a P 5>r 3 IC(i)-I, (27) 

i 

where 

K(i) = 3[S(i) ■ n(i)]n(i) - S(i). (28) 
The components of the axial vector K(i) can be written in the form 

^ 1 (0 = E^, 1 (0, (29) 

I 

where the second-rank, symmetric, traceless tensor 

D qi = 3n g (i)n/(?) - S q in q (i)ni(i) (30) 

is proportional to the normalized spherical function C 2 ^, y?j). 
Then the components of the symmetrized vector K x (i) are 

K}® = c[d 2 ® s 1 ]; = ce ciy^j^^sUi), (31) 

m,m' 
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where the [D 2 (g) S 1 ]* are the components of an irreducible first-rank tensor prod- 
uct of D 2 and S 1 . The constant C can be derived from atomic spectroscopy 
theory [9] and is equal to 10 1 / 2 . 

Using the relations given above and expressing the scalar product of K(i) and 
I in Eq.(27) as a tensor product of their symmetrized forms K x (i) and I 1 , one 
obtains 

H£ = (F^) 1 • Q 1 = £(-ir(i^)J ■ Qi q . (32) 

i 

The q component of the electronic operator has the form 

(F dtp ) 1 q = g aJ2r; 3 K 1 q (t) = 

i 

VWg a ^ Ciy^r-'Cli^Slii), (33) 

m,m' i 

and the q component of the nuclear operator is 

Q- g = ga P i- q - 

The Q l _ q components can be obtained by means of the Wigner-Eckart theorem. 
We have 

(II\Q]\II) = 9<x P ^=j C\ Yz < I\\Q k P > (34) 

= g^iciYi/cjYAM,, 

where Mi = fi is the magnetic dipole nuclear moment. 

By analogy with Eq.(5), the matrix elements < f3IaJFM F \H.^ p \/3IaJ'FM F > 
can be written as 

< (3IaJFM F \U^\(3IaJ'FM F >= 

(-iy+ J '- F W(IIJJ'- F) < aJ\\{Y dip f\\aJ' >< Q 1 > . (35) 

In order to calculate the contribution of this interaction to the hyperfine split- 
ting and to the constant of this splitting A k=1 , we recouple the matrix elements 
< aLSJ\{¥ dip f\aLSJ' > from the LM L SM s JMj to the LM L SM S coupling 
scheme. Since tensor [C 2 (-#, ip) <S> S 1 ] 1 acts over space and spin subsystems the 
9j-symbol technique leads to 



aLSJ\(F dip )\\aLSJ' >= 10V3g a ^{2J + 1)(2J' + 1) 

mm' 

x < aLS^rfCl^^Sl^WLS > . (36) 



< 



x 



f L L 2 ' 

s s i >^c 2 n \Y q 

J J' I mm ' 
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Using Eqs.(61)-(64) of the Appendix we express the matrix elements of the 
electronic factor in the form 

(aLSM s \ Y^rfCl^^Sld^aLSM's) = 

i 

r'j=ri 

< SM s \Sl\SM' s > J r^CK&utpjDsirJdTu (37) 
Finally for the matrix elements of the electronic dipole operator we obtain 
< aLSJ\{F dip )\\aLSJ' >= 15V2g a-^==^/(2J + 1)(2 J' + 1) 




Ellipsis 
°m m' q^MsmM'g 



x J r^C 2 m (A, V i)Ds(r 1 )dr 1 . (38) 



ri 



One can easily see that the functional < aLS J\(F dip )^\aLS J' > coincides, up 
to the constant, with the functional of the electronic multipole operator 
< aLSJ\(F^ 9 ) k q \aLSJ' > (Eq.(25)) for k = 1. 



3.4 Electron-Nucleus Orbital-Spin Operator 

The operator of the interaction between the electronic orbital momentum and 
the magnetic moment of nucleus Eq.(9) can be written down as a product of 
electronic and nuclear operators and has the form 

H£ = (Fl) 1 • Q 1 = E(-1)W; • Q-r (39) 

i 

The q component of the electronic operator is 

(^ = 2a£r- 3 Lj(i) (40) 

i 

while the q component of the nuclear operator is 

Q—q = 9 a pl-q- 

Using the Racah technique we obtain 

< (3IaJFM F \U^\(3IaJ'FM F >= 

(—1) I+J '~ F W(IIJJ'; F) < aJ\\{¥ L y\\aJ' >< fil^W > • (41) 
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The matrix elements (I I\Q q \I I) of the nuclear operator components Q 1 are again 
taken from Eq.(34). 

The matrix elements of the components of electronic operator in the 
LMlSMs coupling scheme are 

< aLSJ\(F L )\\aLSJ' > = 2a(-l) s ~ L - J+k+l ^ (2J + 1)(2J' + 1) 

x W(LLJJ'; kS) < aLS\ ]T rfL\{i)\aLS > . (42) 

i 

And now we may represent the right side matrix elements in Eq.(42)in terms 
of the charge distribution density functional. Thus we obtain 

(aLS | ]Tr- 3 Lj(i) | aLS) = [ rfL\{l)p(KK | nj^dn. (43) 

r'j=ri 

Here p{KK\r\\ r[) = p{K\v\\ r[) is the non-local density functional charge distri- 
bution. 

3.5 Electron-Nucleus Contact Operator 

The operator of the electron-nucleus contact interactions Eq.(lO) can be also 
given as 

H£* = (Fcont) 1 ■ Q 1 = Ei-mFcont)] • Qi q . (44) 

Here the q component of the electronic operator is 

(W! = YflbaE%)Sj(0 (45) 
and the g component of the nuclear operator is 

Q- q = gupl- q - 
By analogy with Eq.(41) we obtain 

< [3IaJFM F \H^ nt \f3laJ'FM F >= 

(-iy+J'- F W(IIJJ';F) < aJ\\(f cont ) l \\aJ' >< PI^W > • (46) 

The matrix elements of the components of electronic operator (F^^) 1 in the 
LM L SM S coupling scheme are 

< aLSJ\(F cont )l\aLSJ' > = ^g a(-l) s - L - J+k+1 ] /(2J + 1)(2 J' + 1) 
x W(LLJJ'; kS) < aLS\Y,$(ri)Sl(i)\aLS > . (47) 
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Using the procedure described in subsections 3.2 and 3.3 and the Appendix 
equations, the electronic matrix elements can be represented in terms of the spin 
distribution density functional 

(aLSM s \Y,Kn)S l MaLSM' s ) = {cJJ^/C^I} 

i 

J 5(r 1 )q(KK'\r 1 ;r' 1 ) 1 dr 1 = (SM s \Sl\SM' s ) 

=ri 

x Js(v 1 )Ds(v 1 )dv 1 ^^-^=C^J M lD s (r)\ r=0 , (48) 

ri 

where D s (r)\ r=0 is the spin distribution density functional at the location of the 
nucleus. 



x 

r', =ri 



4 CONCLUSION 

In this paper we suggest an approach for the calculation of the nuclear mul- 
tipole moments. In this approach the electronic matrix elements entering the 
experimentally observed hyperfme interaction constants are expressed in terms 
of the electron density functionals of the charge or spin distributions. The ma- 
trix elements are presented in a form , suitable for numerical implementation. 
These matrix elements are determined for all the electronic operators entering the 
Hamiltonians of the different descriptions of hyperfme electron-nuclear interac- 
tions. Using the irreducible tensor-operators and the Racah and Fano techniques 
as well as the density functional and density matrix formalism, we can express 
these matrix elements as the products of 3j, 6j, or 9j factors with the space 
part-charge or spin distribution density functional. In principle, the charge or 
spin distribution density functionals can be constructed using every relativistic 
or non-relativistic quantum-mechanical or DFT method. 

The matrix elements of the (F el ) k and (F™ 9 ) fe operators of the Hamiltonians 
Hfc Eq.(l) and the operator (Fl) 1 of Hamiltonian H eri Eq.(7) are the charge 
distribution density functionals while the matrix elements of operator (F™ 9 ) k of 
the Hamiltonians Hfc and two operators (F^p) 1 and (Fcont) 1 of the Hamiltonian 
H en are spin distribution density functionals. It is important to point out that 
the charge distribution density functional does not depend on the level of the 
fine structure, this dependence is governed by the j- symbols. One may notice a 
similar behaviour in the case of spin distribution density functional. But in this 
case the normalized functional remains the same within one given spin multiplet. 
And here, too the fine structure level for which we consider its hyperfme splitting 
within the spin multiplet is determined by 3j, 6j or 9j factors Eq.(25), Eq.(38), 
Eq.(48). Of course, the hyperfme splitting energy depends on j-symbols as well 
as on the corresponding density functional. 
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It is straitforward to point out that such a representation is valid for the 
arbitrary degenerate or non-degenerate states of the electronic system and also 
for arbitrary transitions induced by an external magnetic field. Hence, in the 
case of the free atoms this approach may be successfully applied to all kinds of 
experiments determining the nuclear multipole moments. 

The charge and spin distribution density functional may be constructed by 
means of proper relativistic ab-initio results (for instance, the results of paper 
[19]) or DFT calculations, however of the type which contains more precise rela- 
tivistic corrections [20]. One also may use the density functional of relativistic 
corrections given in papers [15, 16]. Of course our approach gives the possibility 
to make the non-relativistic calculations. The application of this approach in 
non-relativistic calculations has an advantage because of the density functional 
methods can be used. We have shown in subsections 3.2 and 3.3 that the the 
matrix elements Eq.(25) and Eq.(38) of the F 2 ^ and F dip operators written in 
terms of the spin distribution density functionals have the same analytical form 
in all cases. Similar results hold for the operators sum F dip + F L and the opera- 
tors sum F mg = + F^ ig whose matrix elements, however are calculated only 
numerically, or appear in the special case of one electron above the closed shell 
[7]- 

5 APPENDIX 

DENSITY FUNCTIONALS OF CHARGE AND 
SPIN DISTRIBUTIONS 

5.1 Charge Distribution Density Functionals 

The first order density matrix of a iV-electronic system in state K described by 
a wave function ^k(ti, tn), eigenfunction of the operators S 2 and S z , has the 
form [10, 21] 

P {KK\t i] t' 1 )=N J ^ K (r 1 ,T 2 ,...,T N )^* K (r[,T 2 ,...,T N )dT 2 ...dT N , (49) 

where Tj = (rj,<7j), r iy (1 < i < N) being the position vector and <jj the spin 
variable of the i-th electron, and dri = dYidai. 

The corresponding first order density function is defined by the expression 

p(KK\T 1 )=p(KK\T^)k~r i =p(KK\T 1 -,T 1 ). (50) 

The first order transition density matrix between states K and K' has the 
form 

Vk(ti,t 2 , ...,t n )^* k ,(t[,t 2 , ...,T N ),dT 2 ...,dT N . (51) 
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For the one-density transition function we have a similar expression, which 
follows by analogy with Eq.(50). 

After separation of the space and spin variables, the one-density matrix takes 
the form 

P(r 1 ;r[)= £ p^'(r 1 ;r' 1 ) 7 ( (Jl )7 , *K), (52) 

7,7' 

where the p 1 ' 7 ' (r^, r[) are the space components of one-density matrix and the 
7(a) (7 = a, (3) are the electron spin wave-functions. 

The space density matrix of charge distribution has the form 

p(ri;ri)= / pfa; iftdn = p a ' Q (r i; r[) + p^(n; ri) (53) 

r[=n 

and the space density function of charge distribution 

p(n) = p(n; r;) r;=ri = p«( ri ) + /(n)- (54) 

Here p(ri) is local and p(ri; r' x ) non-local electron density functional. 

The matrix elements of a sum of identical one-particle operators, can be writ- 
ten down as 

{K\Y,m\K') = ^K'(r 1 ,...,r N )\^m\^K(r 1 ,...,r N )) 

i i 

= J F^piKK'ln-r'^dn. (55) 
T[=n 

The one particle operators F(i) do not depend of spin variable we have 

(KlYsHm')^ J HlMKK'ln-yjdru (56) 
i r' 1=ri 

and if in addition the F(i) is the operator of multiplication then 

(K\J2F(i)\K') = J F(l)p(KK'\ ri )d ri . (57) 



5.2 Spin Distribution Density Functionals 

In a spin eigenstate (with eigenvalues S(S + 1) and M ) the spin distribution 
non-local density functional, can be written in the following form: [11, 12, 21] 

9(r i; ri) = l/2[p^(r i; r[) - p^(r i; ri)]. (58) 
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The spin distribution matrices q( M \ri, r[) for different spin eigenstates (with 
M — S, S — 1, — S) can be expressed in terms of the normalized spin distribu- 
tion matrices -Ds(ri; r^), which are independent of M 

9 (M) (ri; r[) = q(KK\ ri ; r[) = ™q(K 7f|r i; r'J = MD s (r i; r[), (59) 

with: 

D s {v 1 -T , l ) = ^q(KK\v l -v' l ). (60) 

Here K is the index of the spin state corresponding to (S z ) = M, K is the index 
of the state corresponding to the maximal value of M — S. 

Using the general definition for the transition spin distribution density func- 
tional we can write [10]-[12] 

q(KK'\T 1 -y i ) 1 m = I S^piKK'lr^r^. (61) 

a[=a 1 

Here the S^JJ) (m = 0, ±1) are symmetrized components of the i-electron spin 
operator S(i). Similar formulae can be written for the transition local spin- 
distribution density functional setting ri = r^. 

The spin-distribution functional for transitions between states K(SM S ) and 
K'(S'Ms') satisfy the relations 

q(KK'\ r l5 v\f m = {C^, /C s s A %} q{KK'\ n; (62) 

where C^^^ are Clebsch-Gordon coefficients. By analogy with Eq.(61) one 
obtains 

q{KK'\v 1 -r' l )] n = j S^piKK^r^. (63) 

a[=ai 

Using the Wigner-Eckart theorem we have 

q(KK'\r ^r'Jl = {Cj 8 ^, g / C§ J #} q(KK'\ r l5 r'J 

= KSMslSUSM^Ds^r',), (64) 

where the (m = 0, ±1) are symmetrized components of the total spin S. 
Equation(64) is a generalization of the equation(59) for a transitions between 
different spin states. 
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